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Extending earlier work, we provide predictions obtained with the Monte Carlo generator PROPHECYdf for 
the decays H —+ ZZ/WW 4Z including the complete electroweak 0{a) corrections and some higher-order 
improvements. The gauge-boson resonances are described in the complex-mass scheme. Here, particular attention 
is paid to a comparison of different final states with identical charges, such as e^e~ fi'^fi~ and jj,'^ fi~ fi^ . 



1. INTRODUCTION 

The primary task of the LHC will be the de- 
tection and the study of the Higgs boson. If it is 
heavier than 140 GeV, it decays dominantly into 
gauge-boson pairs, i.e. into 4 fermions. These 
decays offer the largest discovery potential for a 
Higgs boson with a mass Mh ^ 130 GeV [J, and 
the decay H ^ ZZ ^ 4/ will allow for the most 
accurate measurement of Mh above 130 GeV 
At an e"'"e~ linear collider, these decays will en- 
able measurements of the corresponding branch- 
ing ratios and couplings at the few-per-cent level. 

A kinematical reconstruction of the Higgs bo- 
son and of the virtual W and Z bosons requires 
the study of distributions. Thereby, it is impor- 
tant to include radiative corrections, in particular 
real photon radiation. In addition, the verifica- 
tion of the spin and the CP properties of the Higgs 
boson relies on the study of angular and invariant- 
mass distributions j3l4j . As a consequence a 
Monte Carlo generator for H ZZ/WW -> 4/ 
including electroweak corrections is needed. 

In the past, only the electroweak 0{a) correc- 
tions to decays into on-shell gauge bosons H — > 
ZZ/WW [S] and some leading higher-order cor- 
rections were known. However, in the threshold 
region the on-shell approximation becomes unreli- 
able. Below the gauge-boson-pair thresholds only 
the leading order was known until recently. 



PR0PHECY4f is a recently constructed 
Monte Carlo event generator for H ZZ/WW — > 
4/ that includes electroweak corrections as well 
as some higher-order improvements. Since the 
process with off-shell gauge bosons is consistently 
considered without any on-shell approximations, 
the obtained results are valid above, near, and be- 
low the gauge-boson pair thresholds. Parallel to 
our work, another generator for H — > ZZ ^ 4/, 
including electromagnetic corrections only, has 
been introduced in Ref. [J]. 

In this note we briefly describe the structure 
of PR0PHECY4f and of the underlying calcula- 
tions. Moreover, we extend the numerical results 
of Ref. ^ obtained with this code by paying par- 
ticular attention to final states that differ only in 
the generation quantum numbers of the leptons. 

2. CALCULATIONAL DETAILS 

We have calculated the complete electroweak 
0{a) corrections to the processes H ^ 4/. This 
includes both the corrections to the decays H 
ZZ Af and H WW ^ 4/ and their interfer- 
ence. The calculation of the one-loop diagrams 
has been performed in the conventional 't Hooft- 
Feynman gauge and in the background-field for- 
malism using the conventions of Refs. |S| and Pj, 
respectively. The masses of the external fermions 
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have been neglected whenever possible, i.e. every- 
where but in the mass-singular logarithms. 

For the implementation of the finite widths 
of the gauge bosons we use the "complex-mass 
scheme" , which was introduced in Ref. for 
lowest-order calculations and generalized to the 
one-loop level in Ref. In this approach the 

W- and Z-boson masses are consistently consid- 
ered as complex quantities, defined as the lo- 
cations of the propagator poles in the complex 
plane. The scheme fully respects all relations that 
follow from gauge invariance. A brief description 
of this scheme can also be found in Ref. ^21 ■ 

The amplitudes have been generated with Fey- 
nArts, using the two independent versions 1 and 
3, as described in Refs. ^31 and ^1], respectively. 
The algebraic evaluation has been performed in 
two completely independent ways. One calcu- 
lation is based on an in-house program imple- 
mented in Mathematica, the other has been com- 
pleted with the help of FormCalc ^Hj- The ampli- 
tudes are expressed in terms of standard matrix 
elements and coefficients of tensor integrals as de- 
scribed in the appendix of Ref. ^Hl ■ 

The tensor coefficients are evaluated as in 
the calculation of the corrections to e+e~ — > 
4 fermions |llll7j . They are recursively reduced 
to master integrals at the numerical level. The 
scalar master integrals are evaluated for com- 
plex masses using the methods and results of 
Refs. |18ll9l20j . UV divergences are regulated di- 
mensionally and IR divergences with an infinitesi- 
mal photon mass. Tensor and scalar 5-point func- 
tions are directly expressed in terms of 4-point 
integrals Tensor 4-point and 3-point in- 

tegrals are reduced to scalar integrals with the 
Passarino-Veltman algorithm [221 as long as no 
small Gram determinant appears in the reduc- 
tion. If small Gram determinants occur, two al- 
ternative schemes are applied |2S1- One method 
makes use of expansions of the tensor coefficients 
about the limit of vanishing Gram determinants 
and possibly other kinematical determinants. In 
the second, alternative method we evaluate a spe- 
cific tensor coefficient, the integrand of which is 
logarithmic in Feynman parametrization, by nu- 
merical integration. Then the remaining coeffi- 
cients as well as the standard scalar integral are 



algebraically derived from this coefficient. The 
results of the two different codes, based on the 
different methods described above are in good nu- 
merical agreement. 

Since corrections due to the self-interaction 
of the Higgs boson become important for large 
Higgs-boson masses, we have included the domi- 
nant two-loop corrections to the decay H VV 
proportional to G^M^ in the large-Higgs-mass 
limit which were calculated in Ref. ,24^. 

The matrix elements for the real photonic 
corrections are evaluated using the Weyl-van 
der Waerden spinor technique as formulated in 
Ref. r^SI and have been checked against results 
obtained with Madgraph |26j. The soft and 
coUinear singularities are treated both in the 
dipole subtraction method following Ref. p7] 
and in the phase-space slicing method following 
Ref. |2H| • For the calculation of non-coUinear-safe 
observables we use the extension of the subtrac- 
tion method introduced in Ref. [221 ■ Final-state 
radiation beyond 0{a) is included at the leading- 
logarithmic level using the structure functions 
given in Ref. [30] (see also references therein). 

The phase-space integration is performed with 
Monte Carlo techniques. One code employs a 
multi-channel Monte Carlo generator similar to 
the one implemented in RacoonWW and 
Coffer'^^ (29I31| . the second one uses the adap- 
tive integration program VEGAS [221 ■ 

3. NUMERICAL RESULTS 

We use the G^ scheme, i.e. we define the elec- 
tromagnetic couphng by ac^ = v^G^M-^s^/tt. 
Our lowest-order results include the 0(a)- 
corrected width of the gauge bosons. In the dis- 
tributions the photon has been recombined with 
the nearest charged fermion if the invariant mass 
of the photon-fermion pair is below 5 GeV. More 
details about the setup, all input parameters, and 
more detailed results are provided in Ref. j6,. 

In Table Q the partial decay width includ- 
ing 0{a) corrections is shown for different de- 
cay channels and different values of the Higgs- 
boson mass. In parentheses the statistical er- 
ror of the phase-space integration is indicated, 
and 5 = F/Fo — 1 labels the relative corrections. 
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Table 1 

Partial decay widths for H — > 4 leptons including 0{a) and 0{G^^M^) corrections and corresponding 
relative corrections for various decay channels and different Higgs-boson masses (taken from Ref. 6 ). 





MH[GeV] 


140 


170 


200 


H ^ 




r[MeV] 


6[%] 


r[McV] 


6[%] 


r[McV] 


S[%] 


e~e+/x~/x+ 


corrected 
lowest order 


0.0012628(5) 
0.0012349(4) 


2.3 


0.020162(7) 
0.019624(5) 


2.7 


0.8202(2) 
0.78547(8) 


4.4 


1-1+1-1+ 
I = e, n 


corrected 
lowest order 


0.0006692(2) 
0.0006555(2) 


2.1 


0.010346(3) 
0.010074(2) 


2.7 


0.41019(8) 
0.39286(4) 


4.4 


Ve^+jl-Vfj, 


corrected 
lowest order 


0.04807(2) 
0.04638(1) 


3.7 


4.3109(9) 
4.0610(7) 


6.2 


12.499(3) 
11.907(2) 


5.0 


I = C, fl 


corrected 
lowest order 


0.04914(2) 
0.04738(2) 


3.7 


4.344(1) 
4.0926(8) 


6.1 


14.133(3) 
13.458(2) 


5.0 



The first two channels, e-e+^- ^+ and 1-1+1-1+ , 
I = e, fl, result from the decay H ^ ZZ — > 4/. 
The partial widths for only electrons or muons in 
the final state are equal in the limit of vanish- 
ing external fermion masses, since for collinear- 
safe observables, such as the partial widths, the 
fermion-mass logarithms cancel. The width for 
H 1-1+1-1+ is smaller by about a factor 2, be- 
cause it gets a factor 1/4 for identical particles in 
the final state and it proceeds in lowest order via 
two Feynman diagrams that are related by the 
exchange of two outgoing electrons and that have 
only a small interference. The channel Vce+ fi-pfj, 
results from the decay H WW 4/, while the 
last channel 1/11+ 1- vi receives contributions from 
both the decay into W and into Z bosons. 

As can be seen from the table, the relative 
corrections to e-e+fi-fi+ and 1-1+1-1+ practi- 
cally coincide, demonstrating that the corrections 
are not significantly influenced by the interfer- 
ence between the lowest-order diagrams that ex- 
ist for 1-1+1^1+ . Similarly, there is hardly any 
difference between the corrections to VcG+ijl-v^i 
and vil+l-Di. Plots showing the dependence of 
the partial decay widths on Mh in the range 
Mh = 120-700 GeV, together with the corre- 
sponding corrections, can be found in Ref. [H]. 

As examples for differential cross-sections we 
consider distributions in angles between outgoing 
charged leptons. We note that the applied pho- 
ton recombination does not significantly change 
the angular distributions, because recombining 
soft or coUinear photons with fermions does not 



change their directions significantly. 

In the decay H vil+l'-vi^ neither the Higgs- 
boson nor the W-boson momenta can be recon- 
structed from the decay products. Thus, angular 
distributions can only be studied upon including 
the Higgs-production process. If the Higgs bo- 
son was, however, produced without transverse 
momentum, or if the transverse momentum was 
known, the angle between 1+ and I'- in the plane 
perpendicular to the beam axis could be studied 
without knowledge of the production process. We 
define the transverse angle between 1+ and I'- in 
the frame where Ich.t = as 

COS0;+;-.T = J, TTj r, 

|k/,T|m,,pl 

sgn(sin0,+,-_T) = sgnle^ • (k;,T x kj/^x)}, 

where kj(;) ^ are the fermion momenta transverse 
w.r.t. the unit vector e^, which could be iden- 
tified with the beam direction of a Higgs pro- 
duction process. The corresponding distribution, 
together with the influence of the corrections, is 
shown in Figure U] For Mh > 2Mz the interfer- 
ence between the WW- and ZZ-initiated diagrams 
turns out to be negligible, and the corrections 
look similar for the different final states, in partic- 
ular for those that result from virtual W bosons. 

Finally, we investigate the distribution of the 
angle between like-sign leptons in the decays 
H — > e~e+/i^/i+ and H ii^ fi+ fi^ fi+ in the 
rest frame of the Higgs boson. Figure [5] shows 
that the corrections tend to reduce the well- 
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dr 



d0i+i-,T 

0.07 



MeV 



deg 
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Figure 1. Corrected distribution in the transverse angle between the charged leptons (l.h.s.) and cor- 
rections (r.h.s.) in the decays of the type H vil^l'~vi for Mh — 200 GeV. We define "| interference! 
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Figure 2. Corrected distribution in the angle between I and /' in the decays H I l^l' (l.h.s.) 
and corrections (r.h.s.) for A/h = 200 GeV. We define "| interference! = !2(/x~/i"'"/z~/i+) — (e^e^/z^/i"'")!". 



known enhancement in forward direction. In- 
terference terms in the two Born diagrams for 



H ^ /i /i+ are very small for Afn > 2Afz, 

and the relative corrections for the final states 
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/i /i^/i ^-i^d e e^/i /i+ practically coincide. 

4. CONCLUSIONS 

We have presented further results of the genera- 
tor PROPHECYdf which calculates the complete 
electroweak 0{a) radiative corrections to the 
Higgs-boson decays H ZZ/WW 4/. Par- 
ticular attention has been paid to different final 
states with identical charges, such as e+e^/i+/^^ 
and fi'^ fi~ fi~^ fi~ . 
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